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Level statistics is discussed for XXZ spin chains with discrete symmetries for some values 
of the next-nearest-neighbor (NNN) coupling parameter. We show how the level statistics of 
the finite-size systems depends on the NNN coupling and the XXZ anisotropy, which should 
reflect competition among quantum chaos, integrability and finite-size effects. Here discrete 
symmetries play a central role in our analysis. Evaluating the level-spacing distribution, 
the spectral rigidity and the number variance, we confirm the correspondence between non- 
integrability and Wigner behavior in the spectrum. We also show that non-Wigner behavior 
appears due to mixed symmetries and finite-size effects in some nonintegrable cases. 
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1. Introduction 

Statistical properties of energy levels have been studied for various physical systems in 
terms of the random matrix theory (RMT). For quantum systems, the RMT analysis has been 
applied to characterize quantum chaos and to investigating the integrability of a system. For 
quantum spin systems, we adopt a definition of integrability by the Bethe ansatz: an inte- 
grable model is exactly solvable by the Bethe ansatz. Since a pioneering work,^ the following 
conjecture has been widely accepted: If a given Hamiltonian is integrable by the Bethe ansatz, 
the level-spacing distribution should be described by the Poisson distribution: 

Ppoiis) = exp(-s). (1) 

If it is nonintegrable, the level-spacing distribution should be given by the Wigner distribution, 
i.e. the Wigner surmise for the Gaussian orthogonal ensemble (GOE): 

^'Wig(s) = Yexp(^-— J . (2) 
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In principle, the conjecture is valid only for thermodynamically large systems. Furthermore, 
there is no theoretical support for the conjecture for quantum systems. However, we shall show 
that it is practically effective for finite-size quantum systems. In fact, the above conjecture has 
been numerically confirmed for many quantum spin systems such as correlated spin systems-'^"'^ 
and disordered spin systems. ^"^^ In the Anderson model of disordered systems, Ppoi('S) and 
-Pwig(s) characterize the localized and the metallic phases, respectively.^^ 

It is important to study statistical properties of energy levels for XXZ spin chains, which 
are related to various important quantum spin chains as well as classical lattice models in two 
dimensions. Let us consider a spin-| XXZ spin chain on L sites with next-nearest-neighbor 
(NNN) interaction 

L L 

n = Ji J2i^!Sf^, + + Ai5|S-+i) + J2 + S]S]^, + A2S|5f+2), (3) 

3=1 j=l 

where S°' = (1/2)(t" and {a^,ay, a^) are the Pauli matrices; the periodic boundary conditions 
are imposed. The Hamiltonian (3) is nonintegrable when the NNN coupling J2 is nonzero, while 
it is integrable when J2 vanishes. Here we also note that it coincides with the NNN coupled 
Heisenberg chain when Ai = A2 = 1. When J2 vanishes the system becomes the integrable 
XXZ spin chain, which is one of the most important integrable quantum spin chains. Here, 
quantum integrability should lead to Poisson behavior as the characteristic behavior of level 
statistics. When J2 is nonzero, the characteristic behavior of level statistics should be given 
by Wigner behavior. Here we note that the NNN interaction gives rise to frustration among 
nearest neighboring and next-nearest neighboring spins, which should lead to some chaotic 
behavior in the spectrum. In a previous research, however, unexpected behavior of level- 
spacing distributions has been found for NNN coupled XXZ spin chains.^ Robust non- Wigner 
behavior has been seen, although the NNN coupled chains are nonintegrable. The non- Wigner 
behavior of level-spacing distributions appears particularly when total {Stot)= 0, and is 
roughly given by the average of Ppoi(.s) and Pwig(s). Similar non-Wigncr behavior has been 
observed for a circular billiard when the angTilar momentum = 0, and for an interacting 
two-electron system with the Coulomb interaction in a quantum billiard when = 0.^^ Here 
we should note that Wigner behavior has been discussed in ref. 7 for some XXZ spin chains 
in sectors of S^ot 7^ 0. 

In this paper we show that finite-size effects and discrete symmetries are important for 
analyzing the level statistics of XXZ spin chains. We show explicitly how the characteristic 
property of level statistics of the XXZ spin chains depends on the NNN coupling, the XXZ 
anisotropy and the system size. It is nontrivial to confirm the conjecture for finite XXZ spin 
chains which have the integrable line and some higher symmetrical points in the parameter 
space. Numerically we discuss that the characteristic behavior of level statistics is determined 
through competition among quantum chaos, integrability and finite-size effects. We confirm 
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the correspondence between non-integrability and Wigner behavior in the spectrum. We also 
discuss why the unexpected non-Wigner behavior has appeared for S^^^^ = in ref. 6. We 
expUeitly consider two aspects such as mixed symmetry and some finite-size effects. Further- 
more, in various cases we show under what conditions non-Wigner behavior appears due to 
the two aspects. It seems to be rare that such various cases of non-Wigner behavior have been 
completely understood. 

There is another motivation for the present research: another unexpected behavior of 
level-spacing distribution has been found for the integrable XXZ chain in ref. 6. The level- 
spacing distribution P{s) has shown a novel peak at s = for the anisotropy parameter 
Ai = 0.5. The appearance of the peak is consistent with the 5^2 loop algebra symmetry 
which appears only for special values of Ai. Here we do not consider A2 since J2 = 0. Let us 
introduce a parameter q through the relation Ai = {q + l/q)/2. When g' is a root of unity, the 
integrable XXZ Haniiltonian commutes with the 5^2 loop algebra. Here the loop algebra is an 
infinite-dimensional Lie algebra, and the dimensions of some degenerate eigenspaces increase 
exponentially with respect to the system size."'^'^'^^ Thus the degenerate multiplicity of the 
non-Abelian symmetry can be extremely large. Furthermore, the s/2 loop algebra is closely 
related to Onsager's algebra^^ through which the Ising model was originally solved for the first 
time. The XXZ spin chains are thus closely related to the most important families of integrable 
systems through the integrable point, while they are also extended into nonintegrable systems 
quite naturally. We may therefore expect that the RMT analysis of the XXZ spin chains is 
important in discussing level statistics for other quantum systems that have some connection 
to an integrable system. 

The organization of this paper is the following. In § 2, we recall some aspects of numerical 
procedure for level statistics. In particular, we explain desymmetrization of the XXZ Hamil- 
tonian. We also remark that in the paper we consider only the XY-like region where |A| < 1 
for A = Ai and A = A2. In § 3, we show how competition among quantum chaos, quantum 
integrability and finite-size effects appears in the level statistics of the NNN coupled XXZ 
spin chains. We evaluate level-spacing distributions, spectral rigidities, and number variances, 
and confirm the RMT correspondence between non-integrability and Wigner behavior in the 
spectrum. We solve the non-Wigner behavior reported in ref. 6 and show how the behavior of 
level statistics changes due to mixed symmetry and finite-size effects. We also show that the 
characteristic behavior of level statistics does not depend on the energy range. This makes 
remarkable contrast to the level statistics of spinless fermions evaluated in the low energy 
spectrum at 1/3 filling,^° as well as to that of the Anderson model evaluated at edge regions 
of spectrum. In § 4 we discuss level statistics for a special case of the integrable XXZ spin 
chain where it has the s/2 loop algebra symmetry. We observe that there remains many spec- 
tral degeneracies in the sector of S^^^ = even after desymmetrization with respect to spin 
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reversal symmetry. Finally, we give conclusions in § 5. 

2. Numerical Procedure 

Let us discuss desymmctrization of the Hamiltonians of the XXZ spin chains. When per- 
forming calculation on level statistics, one has to separate the Hamiltonian matrices into some 
sectors; in each sector, the cigcnstatcs have the same quantum numbers. The NNN coupled 
XXZ chains are invariant under spin rotation around the z-axis, translation, reflection, and 
spin reversal. Therefore we consider quantum numbers for the total ('S'fot)) the total mo- 
mentum Ktot, the parity, and the spin reversal. However the total momentum Ktot is invariant 
under reflection only when Ktot = or vr. Thus the desymmctrization according to parity is 
performed only when Ktot = or vr. Similarly, Stot is invariant under spin reversal only 
when — 0- Thus the desymmetrization according to spin reversal is performed only when 



It is convenient to use a momentum-based form for the Hamiltonian when we calculate 
eigenvalues of the NNN coupled chains. To obtain the form, we perform the Jordan- Wigner 
and the Fourier transformations on the original spin Hamiltonian. Some details are explained 
in Appendix A. To calculate the eigenvalues, we use standard numerical methods, which are 
contained in the LAPACK library. 

To find universal statistical properties of the Hamiltonians, one has to deal with unfolded 
eigenvalues instead of raw eigenvalues. The unfolding method is detailed in refs. 6 and 12. 

To analyze spectral properties, in this paper, we calculate three quantities: level-spacing 
distribution P{s), spectral rigidity A3(/), and number variance The level-spacing distri- 

bution is the probability function P{s) of nearest-neighbor level-spacing s = Xj+i — Xi, where 
Xj's are unfolded eigenvalues. The level-spacing distribution is calculated over the whole spec- 
trum of unfolded eigenvalues unless we specify the range. The spectral rigidity is given by 



where Nu{e) = T,i9{£ — £j) is the integrated density of unfolded eigenvalues and Og,, denotes 
an average over eq. The average is done on the whole spectrum except about 15 levels on each 
side of the spectrum. The expression of A3(^) gives the least square deviation of N^^s) from 
the best fit straight line in an interval I. The number variance is given by 



where denotes an average over eq-^^ The average is done on the whole spectrum except 
about 10 levels on each side of the spectrum. 

We calculate the spectrum for the 18-site chains with NNN couplings. The matrix size 
is given by the following: 1387 x 1387 for S^^^ = and i^tot = (Here desymmetrization is 







4/19 



J. Phys. Soc. Jpn. Full Paper 

performed except for spin reversal); 1364 x 1364 for S^^^ = and i^tot = 2it/L, where L is 
the number of sites (Here desymmetrization is performed also for spin reversal); 1282 x 1282 
for = 1 and Ktot = 0. 

Numerical calculations are performed for the Xy-like region, |A| < 1, where A is the 
anisotropic parameter Ai or A2. It may be interesting to study for the Ising-like region, 
|A| > 1, because there exist Ising-like magnets. For example, CsCoBrs and CSC0CI3 are 
quasi-lD Ising-like antiferromagnets with A ~ 10.^^ For A ^ 1, however, level statistics is 
not reliable because energy spectra have some large gaps relative to A and the above unfolding 
method is invalid. 

3. Next-Nearest-Neighbor Coupled XXZ Spin Chains 

We now discuss numerically in the section that for the XXZ spin chains the characteristic 
behavior of level statistics is determined through competition among quantum chaos, quantum 
integrability and finite-size effects. 

3. 1 Spin reversal symmetry and Wigner behavior for Sf^^ = 

For a sector of S^^^ = we numerically discuss the characteristic behavior of level statistics 
on the XXZ spin chains. Here we note that spin reversal symmetry has not been considered 
explicitly in previous studies of level statistics for various quantum spin chains. In some sense, 
desymmetrizing the Hamiltonian with respect to spin reversal symmetry has been avoided due 
to some technical difficulty. Level statistics has been discussed only for sectors of S^^^ 7^ 0, 
where there is no need of the desymmetrization with respect to spin reversal symmetry. 

Let us show explicitly such a case that Wigner behavior appears for = if we consider 
spin reversal symmetry. In Fig. 1 we have obtained the numerical results for level statistics 
such as the level-spacing distribution P{s), the spectral rigidity A3(Z) and the number variance 
S^(Z), for the sector of i^tot = ^tt/L and S^^^ = 0. Here we note that in the sector parity 
invariance does not exist and we focus on spin reversal symmetry. The numerical results of 
level statistics shown in Fig. 1 clearly suggest Wigner behavior. The curve of the Wigner 
distribution fits well to the data of the level-spacing distribution P(s) in the main panel. The 
plots of the spectral rigidity are consistent with the curve of Wigner behavior especially 

for small I as shown in Fig. 1. It is also the case with the number variance S^(Z). We see 
small deviations of A3(Z) and for large I because of some finite-size effects. We have 

thus confirmed that Wigner behavior appears also in the sector of S^^^ = for the XXZ spin 
chains with the NNN interaction. 

In the inset of Fig. 1, we have shown the level-spacing distribution for which we do not 
perform the desymmetrization with respect to spin reversal. It shows non- Wigner behavior, 
which is similar to that of some other cases of mixed symmetry as we shall show in Fig. 4. 
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Fig. 1. Level-spacing distribution P{s), spectral rigidity As{l), and number variance of the 

NNN coupled chain for L = 18, J2/J1 0.5, Ai = A2 0.5, S^-^^ = 0, isTtot = ^n/L under 
complete desymmetrization. Broken lines correspond to Poisson behavior, and solid curves Wigner 
behavior. The inset shows P(s) of the same system under incomplete desymmetrization. Namely, 
the desymmetrization with respect to spin reversal symmetry is not performed for the inset. 



The spin reversal operation on the spin variable of the jth site is defined by 

Sf^Sf, S^^-S^. (6) 

Here, S"^ = {Sj ± iSj)/2. Let M denote the number of down-spins in a given sector. The 
value of the total spin operator S^^^ is given by S^^^ = L/2 — M. Some details of spin reversal 
operation are given in Appendices B and C. 
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Fig. 2. The diagram of contour lines of r\ for the NNN coupled chain with L = 18 (in the inset, 
L = 16), in the sector of S^^^ — 1 and i^tot = 0- Roughly speaking, the area above the thick long- 
dashed line, Wigner behavior; the area below the thick solid line, Poisson behavior. The points 
(a) and (b) correspond to Figs. 3(a) and 3(b), respectively. 



3.2 Finite-size effects on the level-spacing distribution 

Let us explicitly discuss finite-size effects appearing in level statistics. They are impor- 
tant in the Poisson-like or non- Wigner behavior observed in level statistics for the completely 
desymmetrized XXZ Hamiltonians. There are two regions in which finite-size effects are promi- 
nent: A region where J2 is close to zero and another region where Ai and A2 are close to 1. 
In the former region quantum integrability appears through finite-size effects, and the charac- 
teristic behavior of level statistics becomes close to Poisson-like behavior. In the latter region, 
Poisson-like behavior appears due to the symmetry enhancement at the point of Ai = A2 = 1, 
where the U{1) symmetry of the XXZ spin chain expands into the spin SU{2) symmetry. 

Let us now discuss how the degree of non- Wigner behavior depends on the anisotropy 
parameters, Ai and A2, and the NNN coupling, J2. For simplicity we set Ai = A2 and 
denote it by A, and we also consider the ratio of J2/ Ji- We express the degree of non- Wigner 
behavior by the following parameter: 

j^\P{s) - Py,,,{s)]ds 

where sq = 0.4729 • • • is the intersection point of i-poi(s) and -Pwig('S)-^'^^ We have r] = when 
P{s) coincides with Pwig{s), and rj = 1 when P{s) coincides with Ppoi(s). The diagram of 
contour lines of rj is shown in Fig. 2. We have calculated them for the area —0.98 < A < —0.02, 
0.02 < A < 0.98, and 0.02 < J2/J1 < 1, where A = Ai = A2. 

The contour lines of rj show that a behavior close to Wigner one appears in a large region, 
while the Poisson-like behavior appears in a narrow region along the line of J2/J1 = and 
that of Ai = A2 = 1. The Poisson-like behavior is dominated by finite-size effects and hence 
should vanish when L ^ 00. This expectation is consistent with the suggestion in ref. 7 that 
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an infinitesimal intcgrability-brcaking term (the NNN term of cq. (3) in this paper) would 
lead to Wigner behavior. In fact, it is seen in the inset of Fig. 2 that the region of Wigner 
behavior shrinks for L = 16. Here we remark that the phase diagrams of the ground state^^'^^ 
are totally different from the diagram of contour lines of r]. It is due to the fact that level 
statistics reflects highly excited states rather than the ground state. 

When Ai = A2 = 0.98, Poisson-like behavior appears in level statistics due to some 
finite-size effects. It will be explicitly shown in Figs. 4 and 5. When Ai = A2 = 1, eq. (3) 
coincides with the Heiscnberg chain, which has the spin SU{2) symmetry. Some degenerate 
energy levels at Ai = A2 = 1 can become nondegenerate when Ai and A2 are not equal to 1. 
The difference among the nondegenerate energy levels should be smaller than the typical level 
spacing when Ai and A2 are close to 1. The typical level spacing, which is of the order of 1/L, 
should become large when the system size L is small. Thus, the Poisson-like behavior should 
practically appear in level statistics. We note that for the Heisenberg chain Wigner behavior 
appears in the level-spacing distribution when we desymmetrize the Hamiltonian with respect 
to the spin SU{2) symmetry.^' ^ 



3. 3 Homogeneity of the characteristic behavior of level statistics throughout the spectrum 

Let us discuss that for the XXZ spin chains the characteristic behavior of level statistics 
does not depend on the energy range of the spectrum. In Fig. 3, we show level-spacing distri- 
butions evaluated at points (a) and (b) shown in the diagram of Fig. 2. They are evaluated 
for three different energy ranges. The distributions shown in Fig. 3(a) give Wigner behavior, 
while the distributions of Fig. 3(b) are close to Poisson behavior. 

Let us explain the three different energy ranges shown in Fig. 3. Histograms show the 
level-spacing distributions evaluated for all levels, while bars show those evaluated only for 
the 1/3 of all levels around the center, and crosses for the 10% of all levels located from each 
of the two spectral edges. 

Quite interestingly, the distributions evaluated for the different energy ranges are quite 
similar to each other. It should be typical of frustrated quantum systems. In frustrated quan- 
tum systems, quantum chaotic behavior appears already in the low energy region near the 
ground state.^^ For non-frustrated cases, however, level statistics shows Poisson-like behav- 
ior in a low region. In ref. 20, for example, level statistics is Poisson-like for the low energy 
spectrum of one-dimensional spinless fermions with the nearest hopping, the nearest- and 
next-nearest-neighbor interactions in a 1/3-filling case. The model is related to our model by 
a transformation, while, the 1/3-filling case is not a frustrated case. Furthermore, there is an- 
other example of Poisson behavior. The level statistics of the Anderson model shows Poisson 
behavior even in the metallic phase, if we evaluate it in the edge regions of energy spectrum. 
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Fig. 3. Level-spacing distribution of the NNN coupled chain for L = 18, Ai = A2 = 0.5, S^^^ — 1, 
and (a) J2/J1 = 0.5; (b) J2/J1 = 0.2. Red histograms are for all levels; green bars, 1/3 of all levels 
around the center; blue crosses, 10% of all levels from each of the two edges. Solid and broken 
lines show the Wigner and Poisson distributions, respectively. 



3.4 Two solutions to unexpected non-Wigner behavior: mixed symmetry and finite-size effects 
Unexpected non-Wigner behavior has been reported in ref. 6 for level-spacing distributions 
of the NNN coupled XXZ chains. Let us discuss the reason why it was observed, considering 
both mixed symmetry and finite-size effects. Here we explain mixed symmetry in the following: 
Suppose that the Hamiltonian of a system is invariant under symmetry operations Tj for 
j = 1,2, which are commuting each other. The eigenstates of the Hamiltonian are 

specified by the set of the eigenvalues of all Tj . If we desymmetrize the spectrum with respect 
to a partial set of symmetry operations, we say that the derived spectrum has mixed symmetry. 
For instance, if we only consider symmetry operations Tj for j = 1,2, ... ,m — 1, then the 
contributions of such eigenstates with different eigenvalues of can be mixed in the derived 
spectrum. 

There are two types of non-Wigner profiles reported in ref. 6 for the nonintegrable systems: 
one is given by almost the numerical average of the Poisson and the Wigner distributions, and 
another one is rather close to the Poisson distribution. The profiles of the first type appear 
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Fig. 4. Level-spacing distribution P{s) of the NNN coupled chain for L = 18, J2/J1 = 1, -fi'tot = 0. 
Broken lines, the Poisson distribution; solid curves, the Wigner distribution. 

in various cases,^ such as the case of A2 = 0.5. The profiles of the second type appear in 
particular for the case of Ai ~ A2 ^ 1. We may call the latter Poisson- like behavior rather 
than simple non- Wigner behavior. Both types of non- Wigner distributions have been observed 
in the subspace of S^^^. = 0, which is the largest sector of the Hamiltonian matrix of eq. (3). 
Here we note that the observations in ref. 6 are obtained particularly for S^^^ = 0, while in 
ref. 7 the Wigner behavior is observed for S^^^ 7^ in the level-spacing distributions of similar 
XXZ chains. 

We give level-spacing distributions in Fig. 4 for the four cases: S^^^ = or 1 and Ai = 
A2 = 0.5 or 0.98. The numerical results suggest that the value of S^^^ should be important 
as well as the anisotropy parameters, Ai and A2, in the observed non- Wigner behavior of 
the level-spacing distributions. When Ai = A2 = 0.5, Wigner behavior appears for S^^^ = 1, 
while the non- Wigner behavior was observed for S^^^ = 0. We have also checked that Wigner 
behavior appears for S^^^ = 2. Furthermore, we have confirmed that such Sfo^-dependence of 
the level-spacing distribution is valid for some values of i^tot- Here we have desymmetrized 
the Hamiltonian according to S^^^, -fCtot) smd the parity when it exists, but not to the spin 
reversal. Here we note that the parity invariance exists only for sectors with Ktot = or tt 
when L is even. 

The non- Wigner behavior observed for the case S^^^ = and Ai = A2 = 0.5 shown in 
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Fig. 5. Spectral rigidity A3(Z) and number variance of the NNN coupled chain for L = 18, 

J2/J1 = 1, -f^tot = 0. Broken lines, Poisson behavior; Solid lines, Wigner behavior. In each of the 
four distributions P{s) vanishes at s = 0: there is no degeneracy among energy levels. 

Fig. 4 is due to mixed symmetry. Let us recall that the system is invariant under spin rotation 
around the z-axis, translation and reflection (parity). For S^^^ = 0, the system is also invariant 
under another operation, spin reversal. Here, the spectrum shown in Fig. 4 has not been 
desymmetrized with respect to the spin reversal symmetry due to some technical aspect. 
However, we expect that Wigner behavior will appear for the case 5tot = and Ai = A2 = 0.5, 
if we further perform desymmetrization with respect to spin reversal symmetry. There are two 
reasons supporting it: First, for i^tot = 27r/L, Wigner behavior for S^^^ = has been shown 
in Fig. 1. Here the system is invariant under spin rotation around the z-axis, translation and 
spin reversal, but not for reflection. Secondly, the behavior of level statistics is independent of 
Ktot for a large number of cases we have investigated. For example, the non- Winger behavior 
for the case S^^^ = and Ai = A2 = 0.5 is very similar to that of the inset of Fig. 1. 

The Poisson-like behavior for the case Ai = A2 = 0.98 is dominated by flnite-size effects. 
In fact, for S^^^ = 1 of Fig. 4, the Poisson-like behavior appears when Ai = A2 = 0.98, while 
Wigner behavior appears when Ai = A2 = 0.5. We have conflrmed that such tendency does 
not depend on the value of i^tot^ we see it not only for Ktot = but also for i^tot 7^ 0. 

The observations of the level-spacing distributions can also be conflrmed by investigating 
spectral rigidity ^^{l) and number variance In Fig. 5, and are shown for 

the four cases corresponding to those of Fig. 4. For S^^^ = 1 and Ai = A2 = 0.5, Wigner 
behavior appears. For S^^^ = and Ai = A2 = 0.5, an intermediate behavior appears, which 
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Fig. 6. Level-spacing distribution P{s) of the integrable chain (J2 = 0) for L = 18, Ai = 0.5, 
S^^^ = 0, Ktot = "iir/L. (a) Desymnictrization with respect to spin reversal is not performed, (b) 
Desymmetrization with respect to spin reversal is performed. The inset is spectral rigidity for each 
case. 



is close to the average between Wigner and Poisson behaviors. For Ai = A2 = 0.98, both for 
^tot — ^tot — 1' Poisson-like behavior appears for A3(i) and S^(Z). 

4. Integrable XXZ Spin Chain in a special case 

Let us discuss level statistics for the special case of the integrable XXZ spin chain that 
has the 5^2 loop algebra symmetry. Here we recall that the XXZ spin chain is integrable when 
the NNN coupling J2 vanishes, and also that the 5/2 loop algebra symmetry exists when qis a 
root of unity. Here the anisotropy Ai is related to q through Ai = (g + q~^)/2. For instance, 
when the parameter q is given by exp(i7r/3), we have A = 0.5. 

The level spacing distribution P{s) and the spectral rigidity ^^{l) of the integrable XXZ 
spin chain are shown for Ai = 0.5 in Fig. 6. In Fig. 6 (a) we do not perform the desymmetriza- 
tion according to spin reversal, while in Fig. 6 (b) we plot the level-spacing distribution P{s) 
and the spectral rigidity A3(/) after we desymmetrize the Hamiltonian with respect to the 
spin reversal operation. 

We observe that there still remain many degeneracies associated with the SI2 loop algebra 



12/19 



J. Phys. Soc. Jpn. Full Paper 

symmetry, even after desymmetrizing the Hamiltonian with respect to spin reversal symmetry. 
The level-spacing distribution P{s) has a small peak at s = in Fig. 6 (b). Furthermore, 
the slopes of A3(/) shown in the insets of Figs. 6 (a) and 6 (b) are larger than that of 
Poisson behavior. However, the numerical result does not necessarily give a counterexample 
to the conjecture of RMT. The level statistics might show Poisson behavior, if wc completely 
desymmetrize the Hamiltonian matrix in terms of the sl2 loop algebra symmetry. 

5. Conclusions 

For the finite spin-^ XXZ spin chains with the NNN interaction, we have evaluated charac- 
teristic quantities of level statistics such as the level-spacing distribution, the spectral rigidity 
and the number variance. Through the numerical results we have obtained the following con- 
jecture: When the symmetry of a finite-size system enhances at some point of the parameter 
space, the characteristic behavior of level statistics should be given by Poisson-like behavior 
near some region close to the point. In particular, we have shown that finite-size effects play an 
important role in the characteristic quantities of level statistics for the XXZ spin chains. Here 
they are integrable for J2 = 0, and their U{1) symmetry extends into SU(2) symmetry at the 
point of Ai = A2 = 1. Furthermore, we have also shown that some unexpected non-Wigner 
behavior appears when an extra symmetry exists and is not considered for desymmetrization, 
such as the case of the spin reversal symmetry in the sector of S^^^ = 0. Here we note that 
in some cases extra symmetries depend on some parameters or some quantum numbers. We 
have thus solved completely the observed non-Wigner behaviors of NNN coupled chain for 
S^,, = in ref. 6. 

Acknowledgement 

The authors would like to thank K. Nakamura and T. Kato for useful discussions. The 
authors also thank the Yukawa Institute for Theoretical Physics at Kyoto University. Discus- 
sions during the YITP workshop YITP-W-03-16 on "Quantum Mechanics and Chaos: From 
Fundamental Problems through Nanosciences" were useful to complete this work. The present 
study was partially supported by the Grant- in- Aid for Encouragement of Young Scientists (A) : 
No. 14702012. It was also partially supported by Hayashi Memorial Foundation for Female 
Natural Scientists. 

Appendix A: Jordan- Wigner and Fourier Transformations 

Let us rewrite eq. (3) by 

"H = Ti-i + 7^2 

1=1 1=1 
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where 7ii is the term containing nearest-neighbor couphngs and is the term containing 
next-nearest-neighbor couplings. We define the Jordan- Wigner transformation by 

exp (-iTT^cjcn cf, (A-2) 



cr, = 



-1 



exp j m^c]cj j q, 
i=i 



(A-3) 



(A-4) 



where and Cj are the creation and annihilation operators of fermions on _7th site. Under the 
Jordan- Wigner transformation, Tii is written by 



Ji 
4 

Ji 



L-l 



1=1 



+ ^[2(4cL+i + 4+iCl) + Ai(44cl4+iCl+i - 24+iCl+i)]. (A-5) 



Here, we consider periodic boundary conditions = 0"^*^): 

CL+i = exp -iTT^ crj"c7t crj^^ = exp 



ivrL ivr v--\ _ 1 , 



(7]^^ exp 



3=1 



mL ivr 



exp iTT^a^. C7+ a^+i = exp — - - ^ 



exp 



ITT 



2 



(^f I 



L 



(A-6) 



where we use e^^af = afe^i'^'^. Since af = ci, <t^ = 4 ' ^tot = ~ M)/2 - M/2 
L/2 — M, where M is the number of fermions, we have 

Therefore, T^i is rewritten as 



(A-7) 



Hi 



L + Y^ +(4c|qc|^iq+i - 4c|q 
1=1 



(A-8) 
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Now, we define the Fourier transformation as 



Here, k takes {2tt/L) x (an integer) for odd M and (27r/L) x (a half-integer) for even M, and 
Q <k < 27r. After the Fourier transformation we have, 

Hi = ^^iL + Ji y^(cos k - Ai)c|.Cfc 

_ JiAi ^ <5fcj+fe2,fc3+fc4exp[-i(A:2-fe4)]4i42^'=3Cfe4> (A-H) 

^1,^2,^3,^4 



where 



{1 when ki + fc2(niod L)=ks + A:4(mod L). 
(A- 12) 
otherwise. 



Considering combination of fc's, we rewrite eq. (A-11) as 
Hi = Ji Ai (^^-M^ +J1YI '^^Cfe 



2JiA 



- XI '^*;i+ik2,fe3+fc4 [cos(A;2 - A;4) - cos(A;2 - A;3)]4i42^'fc3Cfe4- (A-13) 



L 

ki <fe2 
fc3>fc4 



In the same way, 7^2 is rewritten as 



+ X 4i+fc2,fc3+fe4 [cos(A;i + ks) + cos(A;2 + k^) - cos(A;i + ^4) - cos(A;2 + ^3)] cli^kz^kaCki 



L 

^3 >fe4 



2J2A2 



X <5fei+fe2,fe3+fe4 {cos [2 (A;2 - A;4)] - cos [2 {k2 - k^)]} cl^cl^Ck^Ck^. (A- 14) 



L 

ki </c2 
fc3>fc4 

Appendix B: Spin reversal symmetry on momentum-based fermions 

Let us find the momentum-based expression of the mapping corresponding to the spin 
reversal transformation (5^ — Sj, Sj — > —Sj). According to eqs. (A-2), (A-3), (A-9), and 
(A-10), 

1 ^ / i-i \ 

Cfe = ^exp(-iA;Z)exp j -iTr^a-^- cr+ j , (B-1) 



4 = exp(iA;0 exp ( ^/ j " 



(B-2) 
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Under the transformation af^ af , eqs. (B-1) and (B-2) is transformed as follows. 



Cfe 



^ ^ ^ exp(-iA;/) exp |^-i7r ^ aja^ j cT; 

L ' l-l 

1=1 



3=1 

I L ( ^'^ _ \ _ 

= -jf^Z exp[-iA;l - m{l - 1)] exp in V a • C7+ a, , (B-3) 

where Ij is the unit matrix. Now, considering e^'^'' = 1, where Z is an integer, we find that 
exp[— ifcZ — m{l — 1)] = exp[— ifcZ — m{l — 1) + 27ri/] = exp[i(7r — k)l + m] 

= -exp[i(7r- A;)/]. (B-4) 
Therefore, eq. (B-3) is rewritten as 



^ J2 exp[i(7r - k)l] exp ^Itt ^ a^. at j a^ = -cj-fe- (B-5) 



In the same way, we have 

4 -Cn-k- (B-6) 

It is sometimes convenient to use the following in stead of eqs. (B-5) and (B-6): 

4 ^ c^r-fe, Cfc ^ cl_,^. (B-7) 

The Hamiltonians are invariant not only for eqs. (B-5) and (B-6) but also for eq (B-7). The 
form (B-7) has an advantage that we do not need to consider the phase factor —1 that appears 
in eqs. (B-5) and (B-6). 

Making use of the spin reversal operation expressed in terms of the fermion basis (B-7) 
we have desymmetrized the Hamiltonian matrix in the sector S^^^ = with respect to spin 
reversal symmetry. For a given vector with S^^^ = we calculate how it transforms under the 
operation (B-7). If it is not a singlet and transforms into a different vector, then we combine 
the pair into an eigenvector of the operation (B-7). 

Similarly, when we desymmetrize the Hamiltonian with respect to parity, we use the parity 
operation expressed in terms of the fermion basis: 

4 ^ eifegt_^(_i)M-l+L/2-5f„, 2^ ^ ^_ife._^^_^^M+L/2-5f„,_ (g.g) 

Here we recall that M denotes the number of down-spins in the sector of S^„^ = L/2 — M. 
Through the Jordan- Wigner transformation we easily derive the parity operation (B-8) from 
that on the spin variables: — >■ crj_,_i_; , erf — cy^+i-i for Z = 1, 2, . . . , L. 



16/19 



J. Phys. Soc. Jpn. Full Paper 

Appendix C: Spin reversal operation on the vacuum state 

When dcsymmctrizing the Hamiltonian with respect to spin reversal symmetry, it is useful 
to know how the vacuum state transforms under the spin reversal operation expressed in terms 
of the momentum-based fermion operators. Let us denote by |0) the vacuum state where there 
is no down-spin. Under the spin reversal operation it transforms up to a phase factor yl^;, as 
follows 

|0)^Aict^4---ctjO) (C-1) 

Here qj denotes momentum {2Tr/L)j for j = 1,2, . . . ,L, when M is odd, and {2iT/L){j — 1/2) 
for j = 1,2, . . . ,L, when M is even. The phase factor Al is given by 

1 ^ 
= 7172 Yl epexp(-i^jfcp,) . (C-2) 

PG5i j=l 

Here Sl denotes the set of permutations on L elements, ep the sign of permutation P. Fur- 
thermore, we can calculate the phase factor Al as follows 

f (-1)^ for M = 2^ , , 

Al= < (C-3 
\ (-1)^+^ for M = 2^ + 1 

The derivation is given in the following. Let us introduce the following matrix: 

Us = lla^. (C.4) 
We may express the spin reversal operation (6) as follows 

Usclu-^ = -c^_jk, UsCkU-^ = -cl_^. (C-5) 



We thus have 

L 

Us |0) = n |0) = ^1^2 ■■■^l |0) (C-6) 

Applying the Jordan- Wigner transformation and substituting cjs with c|.s through (A- 10), we 
have 

1 ^...^e-Ki^^-^^^^^---+-'=^)<4---4l^ 




^pMi^---^ JO) (C-7) 



Thus we have the expression (C-2) of the phase factor Al. Here we recall that qj = (2iT/L)j 
for j = 1, 2, . . . , L for odd M, and qj = (27r/L)(j - 1/2) for j = 1, 2, . . . , L for even M. We 
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also recall that M denotes the number of down-spins in the sector. 

We now calculate the expression (C-3) for the phase factor A^. We take a vector \v) with 
^tot = as follows. When M is odd, we introduce i = {M — l)/2 and we define \v) by 

b) = (44---4)-(clr--cl2cli)-4io) 

Here Cj denotes with k = (27r/L)j. When M is even, we take £ = M/2, and we define \v) 

by 

1^) ^ (4/24/2 ■ ■ ■ 4-1/2) ■ (^-(€-1/2) ■ ■ ■ ^-3/2^-1/2) 

Here 4+1/2 denotes c|. with k = (27r/L)(j + 1/2). Through the operation (B-5) we show that 
U,\v) = {-1Y+'^Al\v) for M odd and Us\v) = {-1YAl\v) for M even. Thus, we have at least 
U,\v)=±\v). 

Let us show that Us\v) = +\v) for both odd M and even M cases. First we note 
that J/scjcl • • • 4/2l^) ~ '^L/2+i'^L/2+2 ■ ■ ■ '^i/l^)- Second, expanding the vector \v) in terms 
of 4i4o'''4 |0) with 1 < ji < j2 < ■■■ < Jl/2 < L. we show that the coefficient of 
44 ■ ■ ■ ^L/2l^) expansion is equal to that of c|^/2+i'^L/2+2 ' ' ' Therefore we have 

Vs\v) = +\v). 

Thus, we obtain the expression (C-3) for the phase factor Al. 
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